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Dynamics of magnetic flux tubes in close binary stars 
I. Equilibrium and stability properties 
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Abstract. Surface reconstructions of active close binary stars based on photometric and spectroscopic observations reveal non- 
uniform starspot distributions, which indicate the existence of preferred spot longitudes (with respect to the companion star). 
We consider the equilibrium and linear stability of toroidal magnetic flux tubes in close binaries to examine whether tidal effects 
are capable to initiate the formation of rising flux loops at preferred longitudes near the bottom of the stellar convection zone. 
The tidal force and the deviation of the stellar structure from spherical symmetry are treated in lowest-order perturbation theory 
assuming synchronised close binaries with orbital periods of a few days. The frequency, growth time, and spatial structure of 
linear eigenmodes are determined by a stability analysis. We find that, despite their small magnitude, tidal effects can lead to 
a considerable longitudinal asymmetry in the formation probability of flux loops, since the breaking of the axial symmetry 
due to the presence of the companion star is reinforced by the sensitive dependence of the stability properties on the stellar 
stratification and by resonance effects. The orientation of preferred longitudes of loop formation depends on the equilibrium 
configuration and the wave number of the dominating eigenmode. The change of the growth times of unstable modes with 
respect to the case of a single star is very small. 
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1. Introduction 

Owing to their rapid rotation and deep convection zone, 
evolved (sub-)giants and cool main sequence components 
in synchronised close binaries like RS CVn and BY Dra 
systems typically show high levels of magnetic activity 
dStra ssmeier et alJll993l) . e.g., in the form of large starspots 
and enhanced chromospheric and coronal emission in the UV 
and X-ray spectral ranges. Analyses of light curves with in- 
version techniques allow to derive the characteristic proper- 
ties of spots like thei r number, posit i on, size, and temperature 
llEaton & Halll 1 19791: iRodono et alJ Il986l: iBudding & ZeiiiS 
Il987l) . Utilising rotation-modulated, asymmetric spectral line 
profiles of spectro- and spectropolarimetric o bservations, in- 
version techniques like the Doppler Imaging dVoet & Penro j 
1983; Vogt et al. 1987) and Zeeman-Doppler Imaging (Semel' 
1989: Don ati et al. 1992. 1999) reveal information about tem- 



Surface reconstructions of close, fast-rotating binaries 
show extende d starspots which cover large fraction s of the 
stellar surface jO'Neal et alJll998HLanza et alJl200ll) and fre- 
quently occur, in contrast to the case of the Sun, also at inter- 
mediate, high and even polar latitudes. Non-uniform spot dis- 
tributions in longitude indicate the existence of preferred longi- 
tudes (PL), where spots occur more frequent or last longer than 
at other longitudes cHenrv et al. 1995; Jetsu 1996). Preferred 
longitudes are often found to be separated by about 180°, i.e., 
in opposite active quadrants of the hemisphere, but their ori- 
entation with respect to the companion star is not unique. In 
short-period systems with rotation periods < 1 d active regions 



tend to prefer the quadrature longitudes (Zeilik et al 


i990a b. 


1994; 


Olah et al. 


1994; Heckertetal. 1998), 



perature and magnetic field inhomogeneities on the stellar sur- 
face, which are interpreted as cool spots caused by emerg- 
ing magnetic fields. Recent observations are complemented by 
photometric data which cover up to several dec ades and thus 
enable the determination of cyclic variations jRodono et alJ 
[1995, .2003). 



I.e., the 

longitudes perpendicular to the line connecting both stellar cen- 
tres. In systems with longer periods the situation is less clear: 
several systems show PL at the substellar point, its antipode 
or ot her fixed longitudes with respect to the companion star 
(e.g.. lLanza et a l.''l998'.'200l':'Olah et al.''2002'). whereas other 
stars reveal a longitudinal migration of the spotted regions 
with respect to the cor npanion stan which is usually ascribed 
to diff'erential rotation (^Strassmeieil ll994lRodon6 et alJI 19951 
Considering the dependence of PL on the orbital period, 
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iHeckert & Ordwavl ( 1 19951) find a transition from fixed PL for 
orbital periods less that one day to migrating and eventually 
no PL at all for more separated systems with periods exceed- 
ing about 10 days, leading to the suggestion that tidal effects 
might influence the spot distribution. However, observations of 
close binaries with giant components show that spot clusters at 
migrating PL may also occur in systems with longer orbital pe- 
riods (|£erdvu2ina & Tuominen 1998; Berdvu2ina et al. 1999), 
giving rise to the assumption that the formation and orientation 
of clusters are not exclusively dependent on the binary separa- 
tion alone. 

In this and a subsequent paper we study the influence of 
tidal effects on the dynamics and evolution of magnetic flux 
tubes in active binary components, assuming - in analogy to 
the case of the Sun - that starspots are formed by erupting flux 
tubes which originate from the bottom of the stellar convec- 
tion zone. The present work delves into the equilibrium and 
linear stability properties of flux tubes stored in the convective 
overshoot region below the convection zone proper in order to 
examine whether the influence of the companion star is able to 
trigger rising loops at preferred longitudes which may penetrate 
to the convection zone above. The decomposition of small per- 
turbations of an equilibrium flux tube in terms of eigenmodes 
and the determination of the corresponding eigenfrequencies 
are accomplished b y a linea r stability analys is which ext ends 
earher treatments bvlFerriz- Mas & Schiissle^ 

ill 993[ll 99.4 and 

iHolzwarth & Schiissled (120 00). In contrast to rotating single 
stars, in binary stars the existence of the companion star breaks 
the rotational symmetry. Consequently, the dynamics of flux 
tubes is not only subject to the tidal force but is additionally 
affected by the deviation of the stellar structure from spherical 
symmetry. The linear stability analysis deals with the response 
of equilibrium flux tubes to small perturbations inside the over- 
shoot region, whereas the growing displacements of unstable 
loops have to be followed by non-linear numerical simulations. 
The evolution of flux loops rising through the convection zone 
and the resulting surface distribution upon their eruption is de- 
scribed in a following paper. 

In Sect. 121 we briefly summarise the underlying (solar) ac- 
tivity model based on erupting magnetic flux tubes and outline 
the binary model. Sect.|3]describes the equilibrium properties 
and Sect. 0] gives the results of the linear stability analysis of 
toroidal flux tubes. Sect.|5]contains a discussion of the results 
and Sect. Ogives our conclusions. 



penetrate from the convection zone above. While Hale's po- 
larity rules and the East- West-orientation of bipolar sunspot 
groups suggest equilibrium flux rings parallel to the equato- 
rial plane, overshooting gas motions cause perturbations which 
can tri gger the onset of an undulato r y (Parker-type) instabil- 
ity ( Spruit & van Ballegooiienlll982t iFerriz-Mas & Schiissled 
11993.. 1995r If a tube segment is lifted upward, a net mass 
downflow from the crest increases the density contrast with re- 
spect to the environment. Once a critical magnetic field strength 
is exceeded the resulting buoyancy surpasses the opposing 
magnetic curvature force, the loop continues to rise through 
the convection zone , and eventual ly erupts at the stellar sur- 
face dParker... 1 955.: .Moreno-Insertis^ 1 98 6). Upon emergence it 
causes the various signatures of magnetic activity like active re- 
gions, dark spot groups, bright faculae, chromospheric activity. 
X-ray bright coronal loops, and flares. The buoyancy-driven 
mechanism implies a crucial dependence on the stellar strat- 
ification, viz. the superadiabaticity, 6 - V - Vad. In the case 
of the Sun, magnetic flux tubes with field strengths of about 
10^ G are required to obtain theoretical results which are in 
accordance with observational constrains like flux emergence 



latitudes, tilt angles, and asymmetries of bipolar st 


)ot groups 


jMoreno-Insertisll 1 9921: ID' Si 


va & Choudhurill 19931 


iFan et alJ 


'1994: 'Schussler et alJll994l 


Caligari et alJll995l Il998l). The 



quantitative elaboration of the applied flux tube model is also 
consistent with t he emergence of pol ar spots on rapidly rotat- 
ing young stars dOranzer et a l.' 2000) and with the conspicu- 
ous decline of coronal X-ray emission in giant stars a cross the 
'coronal dividing line' ( IHolzwarth & SchussleJ200ll) . 



We consider the magnetic flux in the convection zone in the 
form of isolated magnetic flux tubes, which are embedded in a 
field-free environment. The stellar convection zone is treated 
as an ideal plasma with vanishing viscosity and infinite con- 
ductivity, implying the concept of frozen-in magnetic flux. All 
calculations are carri ed out in the frarnework of the 'thin- flux- 
tube approximation' jRoberts & Webblll978l:ISpruill98ll) : the 
diameter of a circular flux tube is assumed to be small with re- 
spect to the other relevant length scales of the problem like the 
pressure scale height, radius of curvature, or the wave length of 
a perturbation along the tube. As a consequence of very short 
signal travel times it remains in lateral (total) pressure equilib- 
rium with its environment. 



2. Model assumptions 

2.1. The 'flux tube paradigm' 

It is widely accepted that sunspots and bipolar spot groups 
are caused by erupting magnetic flux tubes which origi- 
nate from the base of the solar convection zone. The mag- 
netic field is believed to be amplified in the rotational shear 
layer (tachocline) near the base of the convection zone and 
stored in the for m of toroidal flux tubes in the convective 
overs hoot region jMoreno-Insertis et"al]ll992t ISchiissler et alJ 
I1994I) . This stably stratified subadiabatic layer at the inter- 
face to the radiative core is generated by gas motions which 



2.2. Binary model 

Observations indicate that the orbital motion of binary sys- 
tems with periods T < lOd is effectively synchronised with 
stellar rotation and typically exhibit very small ecc entricity 
dPuauennov & Mayorlll99ll: IStrassmeier et alJll993l) . Hence, 
we consider a close, detached, synchronised and circularised 
binary system, for which the orbital and stellar rotation periods 
are equal and the spin axes of the components are perpendicu- 
lar to the orbital plane. The active star with mass and the 
companion star with mass M^o - qM-^ move on circular orbits 
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with a period of a few days. According to Kepler's third law, 
they are separated by the constant distance 

,1/3 .7^.2/3 
(d) ■ 



\Ro) \Mq 



Our reference binary system consists of two solar-type stars 
{Mi, - Mco - 1 Mq, i.e. stellar mass ratio q - \) with a period 
of r = 2d. The active star is described by a perturbed sin- 
gle star (here, standard solar) model, whereas the companion 
star is considered as an idealised point mass. In the case of the 
Sun, the differential rot ation profile affects the stabilit y proper- 
ties only slightly (e.g., Ferriz-Mas & SchiissleJ[l993b . Owing 
to tidal interactions and synchronisation, the diff'erential rota- 
tion in close binary stars is much we aker than in the solar case 
( lDonahuell993tlRudiger et alJl998h . so that we assume solid- 
body rotation. 

For systems with orbital periods of a few days, tidal effects 
are sufficiently small to be treated in lowest-order perturbation 
theory. Using spherical coordinates in a co-rotating frame of 
reference (radius r, latitude A, azimuth 0, with the companion 
star in the direction (p - 0), the effective potential, ^PefF = *P* + 
^Ftide, is the sum of the spherical gravitational potential of the 



star, *I'*(r), and the potential of the tidal perturbation 



'V^i,(r,(l>,A) 



2r 



-q cos A cos ; 



cos A - q 



(2) 



where e = r/a is the expansion parameter and G the gravita- 
tional constant. The potential *Ptide gives rise to the tidal accel- 
eration 

gtideir, (j). A) = e'g-kir) [e^ + 3q (Cr ■ Ca) Ca 

-(I + q) (Cr ■ en)en] + 0{6') , (3) 

where Cr, en, and Ca are unit vectors according to the directions 
indicated in Fig. ^ Following Eq. 0, the component of the 
tidal acceleration in azimuthal direction e^, for example, is 



(glide -e^) = --gi,{r)e^qcosAsm2(l) + 0(^£*) . 



(4) 



Our model includes the deviation of the stellar structure 
from spherical symmetry arising from the tidal acceleration 
g-fide- For a not to strongly deformed star, we suppose that the 
stellar structure is still described appropriately by functions of 
the effective potential, ^eff, only, i.e. f(r, (p. A) = /(*Peff), where 
/ stands for any stellar quantity (like gas pressure, density,. . .). 
This permits an analytical expression for the tidal deformation, 
namely 



fir,<p,A) = /CP*+*Pt,de) 

Hf GMi, 



d>P 



= /o 1 + 



= fn 1 1 + 1- r — cos 20 



(5) 



' Equations J2j and J3j also contain the centrifugal potential and 
acceleration, respectively, although these contributions do not strictly 
belong to tidal interactions but occur in rotating single stars as well. 




Fig. 1. Geometry of the binary system. The origin of the co- 
rotating frame of reference lies in the centre of mass, CM. The 
stability analysis is carried out in spherical coordinates (r, (p. A) 
with respect to the centre of the active star, . 



with the dimensionless coefficients 



1 -I- -^|cos' 



A — q 



-qcos 



(6) 



(7) 



and the scale height Hf - -(d ln//dr)"' . Equation (|5} includes 
the rotational flattening and the tidal (ellipsoidal) deformation, 
with tidal bulges centred on the line connecting the two stars. 
The deviation of a quantity. A/ = f{r, <p, A) - /o(r), from the 
single star stratification, /o, depends on its local scale height, 
Hf. Both values, /o and Hf, are taken from the spherical model 
at radius r. The quality of the approximation can be estimated 
by considering the ratio between the geometrical deformation, 

6r 



- r + rcos2(f> + 0(e ) , 



(8) 



and the smallest scale height Hf of all quantities /. Since in 
the overshoot region, the supposed storage location of mag- 
netic flux, the scale height of the superadiabaticity, Hg, is con- 
siderably smaller than of any other quantity (see Tab. [0, the 
tidal influence on the stability p roperties of flux tubes is mainly 
controlled by the variation of 6 ( iHolzwarth & Schussleii200(]|) . 

For our reference system (T = 2 d, M^, - Mq and q = 1), 
Eq. Q yields a separation a ~ 87?0 and < 10"^. The tidal 
effects are of small magnitude only, but they imply a significant 
qualitative difference, since both Eqs. and exhibit a pe- 
riodic azimuthal dependence, which breaks the axial symmetry 
present in the single star problem. 

3. Stationary equilibrium 

We consider magnetic flux rings in mechanical stationary equi- 
librium, which are embedded in the overshoot region and par- 
allel to the equatorial plane. The tubes are non-buoyant and in 
lateral pressure equilibrium with the environment. The inward 
dkected magnetic curvature force is compensated by the out- 
ward directed Coriolis force, which arises from a prograde flow 
inside the tube with a velocity v in excess of the stellar solid- 
body rotation. The equilibrium tube forms a torus with constant 
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toroidal 
flux tube 



deformed 
overshoot region 



Fig. 2. Sketch of a toroidal equilibrium flux tube embedded 
in the tidally deformed overshoot region. The flux ring cuts 
through different equipotential surfaces, which implies a pe- 
riodic variation of the environment along the tube: the values 
of the stratification at the points on the line of centres, Vi 2, 
correspond to deeper layers than those at the quadrature points, 

Gl,2. 



radius of curvature, Rq - cos /Iq, where ro and Aq are the ra- 
dial and latitudinal equilibrium positions, respectively. Since 
the flux ring is situated in a tidally deformed star (Fig. the 
internal density, p{- p^, where pe is the density of the environ- 
ment), depends, according to the approximation in Eq. Q, on 
the azimuth as 



= 1 + ro— + ro— cos Itp + O(e^) . 

PO HpO npQ 



(9) 



This entails a periodic variation of the magnetic field strength. 



and internal flow velocity. 



(10) 



(11) 



where Mao - vo/cao is the Alfvenic Much number, i.e., the 
flow velocity in units of the Alfven velocity, Cao - BqI ^j4npo. 
Quantities with index '0' refer to a comparable axisymmetric 
equilibrium flux tube in a slowly rotating single star, where the 
tidal and centrifugal acceleration as well as the deformation of 
the stellar structure are negligible (e.g., More no-Insertis et alJ 
[1992). For flux tubes in mechanical equihbrium, these quanti- 
ties obey the relation 



Vq - RqQ. 



1 + 



Roil 



> . 



(12) 



Table 1. Parameters at equilibrium depth ro of the reference 
configuration 



equilibrium depth 


' u 


5.07 ■ 10'" cm = Q.13Rq 


gas pressure 


PcO 


= 4.31 ■ 10"dyn/cm- 


pressure scale height 


H pQ 


= 5.52 ■ 10' cm 


density 


PeO 


= 0.15g/cm^^ 


density scale height 


^pO 


= 9.21 • 10" cm 


superadiabaticity 


So- 


: -9.77 ■ 10-^ 


superad. scale height 


Hso 


= 4.43 ■ 10** cm 


rotation 


T = 


2d 


binary separation 


a = 


8.41 Rq 


expansion parameter 


= 


6.53 • lO-'* 


deformation parameters 


r = 


8.16- 10-'' 


(at equator, A = 0) 


f = 


4.90 • lO^" 



In the following, we consider flux rings in the middle of the 
overshoot region. Table^gives values of the characteristic pa- 
rameters of this reference configuration for a model of a 1 Mq 
star. 

4. Linear stability analysis 

4.1. Determination and interpretation of eigenmodes 

Once all quantities of the external stratification are determined 
by the equilibrium depth ro, the stability properties of the flux 
ring depend only on its latitude, Aq, and magnetic field strength. 
Bo- Lagrangian displacements are decomposed in eigenmodes, 

E(0,f) =f(0)exp(/wf) , (13) 
characterised by the eigenfrequency co and the eigenfunction 
f(0) = (^t,^n,^bV (in the co-moving trihedron with the local 
tangent, normal, and binormal unit vectors, [t,n,b], respec- 
tively). For unstable eigenmodes with increasing amplitudes 
T - l/|3(w)| provides the characteristic growth time. The dis- 
placement vector ^ is determined by the hnearised equations of 
motion, 

f " + [M^ + iojM^t) + [Mi; + iuM, - oj^M,) ^ = , (14) 

where primes denote derivatives with respect to the azimuthal 
coordinate, (see AppendixIXt. Taking the tidal force, Eq. Q, 
and the azimuthal variation of the equilibrium quantities, 
Eqs. (|9|l - (II It - into account, the coefficient matrices depend 
periodically on 0: 

M{(p) = M{4> + 7i) 

= Mo + e^McCos20-i-e^MsSin20-HC)(e'^). (15) 

Closed flux tubes require ^{(p + 2n) = ^((p) and ^'((f> + 2n) = 
^'(0), suggesting the ansatz 

f(0) = |(,^)exp(im<^) (16) 



^ |/texp(il0) 



exp (inKj. 



(17) 



In analogy to the single-star problem, a phase factor with the 
azimuthal wave number m has been separated. Due to the <p- 
dependent contributions in Eq. il5\ . the substitution of Eq. J17> 
into Eq. J14> yields the 3-term recursion formula 



■Cm+k^k-l + Cm+k^k + 'R-m+k^k+2 







(18) 
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k=-2 



m=2 
k=0 



k=2 



k=-2 



k=0 
m=l 



k=2 



Fig. 3. An eigenmode in the binary problem is represented 
by a spectrum of coupled wave modes with wave numbers 
n - m ± k,k = 0,2,4, ■■ ■ Due to the weak coupling (oc e^) the 
eigenmode can still be characterised by the wave number m of 
its dominating constituent. The contributions of wave modes 
with k lead to non-axisymmetric modifications, whereas 
wave modes with \k\ > 4 (indicated by dots) are negligible. 



In contrast to the single-star problem, where an eigenmode cor- 
responds to a single azimuthal wave number m, i.e., - const, 
and - for k Q,an eigenmode here consists of a spectrum 
of coupled wave modes n = m+k with ^ = 0, +2, +4, . . ., whose 
amplitudes are recursively related to each other through 
Eq. (I18> . This relation is illustrated in Fig. |3] Since the cou- 
pling is weak [oc 0(e^)], the amplitudes decrease quite fast 
with increasing \k\, so that more distant wave modes (here, with 
1^1 > 4) do practically not contribute to the eigenmode. In most 
cases the amplitude is much larger than ^+2, so that the re- 
sulting eigenmode can still be characterised by the wave num- 
ber m. 

The linear eigenfunctions characterise the structure of 
small displacements of flux tubes in the overshoot region prior 
to their rapid, non-linear rise through the convection zone. 
Owing to the stellar stratification, the most important aspect 
of an eigenfunction is its component in the radial direction, 

- ii ■ Sr): the larger the radial displacements, the larger 
are the changes in the environment of the corresponding tube 
segments, which become located in either more stable or more 
unstable layers. Assuming that larger radial displacements to- 
ward more unstable layers favour the penetration of a growing 
flux loop to the convection zone, we use the envelope, ||, (0)|, 
of the radial component of the (normalised) eigenfunction as 
a measure for the onset of rapidly rising flux loops at azimuth 
(p. Since for eigenmodes with 3(w) < the magnitude of the 
radial envelope increases exponentially in time. 



|H,(,^,f)l = l|.(<^)|exp[-3(w)f] 



(19) 



an unstable flux tube in the overshoot region will eventually 
enter the superadiabatical region. Its actual displacement de- 
pends, however, also on the azimuthal wave number m, the fre- 
quency !R(w), and the phase relation 



ai-g S,(0, t) = m0 + arg|,(0) + ^(w)f 



(20) 



corresponding to an azimuthal propagation of the tube's crest 
along the envelope, Eq. M9\ . with the phase velocity 



Vp(<^) = 1 + - — arg^,(0) 

m \ m af 



(21) 



'< 40 




Fig. 4. Relative deviation of the growth time, (xbin - TsiiO/xsin, 
in parameter domains which have been unstable in the rotating 
single-star problem, i.e., Tsin + 0. 



The radial envelope thus represents a kind of 'statistical 
weight', which is applicable either to a large ensemble of 
comparable flux tubes with similar equilibrium properties and 
phase shifts evenly distributed in the range < arg 3^ < 2n, or 
to an individual flux tube with a growth time much longer than 
the wave frequency, t » Inl^icS). If an eigenmode is dom- 
inated by the wave mode m and only slightly modified by the 
adjacent wave modes m + 2, i.e., ^,-.±2 |/-,o := 1, the radial 
envelope is approximately given by 

\U<P)\ = ||,o+|,-y'^+|,-,-2e-''2^ + ...| 

^ 1 + ||r,2| COS (2<^ + ai) + \^r,-l\ COS (20 - Q'_2) 



(22) 



with a ±2 - argfr,±2- The eigenfunctions thus inherit the prop- 
erties of the underlying problem and exhibit a 7r-periodicity in 
azimuthal direction. Assuming that the shape of \^r{(p)\ repre- 
sents a measure for the onset of rising loops, Eq. M2\ implies 
the existence of preferred longitudes at the maxima of \^rV 

4.2. Growth times and eigenfrequencies 

The influence of the companion star modifies the stability prop- 
erties of flux tubes by breaking the axial symmetry. The result- 
ing quantitative diff'erence in growth times r with respect to a 
rotating single star is very small; as illustrated in Fig.|4] the rel- 
ative deviations in all cases are considerably smaller than 1%. 
Qualitative differences occur, however, for equilibrium config- 
urations which had been stable in the single-star problem. In 
previously stable parameter domains new types of instabilities 
with very long growth times lead to an 'instability background' 
consisting of extended 'plateaus' which are pervaded by nar- 
row 'ridges' (see Fig. |5}. 

The 'plateaus' are due to the tidal force and the azimuthal 
variation of the density along the equilibrium flux ring. Using 
the lateral (total) pressure equilibrium and the variation of the 
magnetic field strength, Eq. ilO\ . it can be shown that the in- 
ternal specific entropy, S, is also a function of (p and varies 
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Fig. 5. Growth times of the 'instability background' due to 
tidal effects. In addition to the Parker-type instabilities there 
are localised ridges with t > 10^ ... 10^ d traversing extended 
plateaus with r > 10^ d. Because of the limited resolution of 
the parameter grid in the figure the very narrow ridge at A ^ 1 ° 
appears to be discontinuous. 

approximately as 

^dS_ Idlnpe 

Cp dcf, ~ /3 dcf> ~ p Hp,, ' 

where yS = ^np/B^ is the ratio between the gas pressure and the 
magnetic pressure, and Hp Q is the local pressure scale height. 
Although we have r/Hpe ^ 9, the azimuthal variation is nev- 
ertheless small compared to the radial variation of the stellar 
stratification, since the relevant field strengths considered here 
lead to values in the regime (3 ~ lO'* . . . 10^ » 1. Due to its 
velocity excess, v, a gas element inside the tube would have 
to change its entropy to preserve the equilibrium conditions. 
But since there is no mechanism available to accomplish this 
entropy change, the azimuthal variation leads to a loss of the 
mechanical equilibrium and to very slow monotonic displace- 
ments of the flux tube out of its toroidal configuration. For fast- 
rotating stars, the internal flow velocity v is small and thus the 
growth times of plateau instabilities are typically several 10^ 
days (up to t ~ 10' d). The effects of the azimuthal density 
variation and tidal force along the tube are negligible because 
the stability properties are dominated by the azimuthal varia- 
tion of the superadiabaticity, 

^ 6(ro)^l + , (24) 

owing to its relative small scale height in the overshoot region 
(with Hg/Hp ^ 0.08, see Tab. [Q. If all azimuthal variations 
other than that of 6 are omitted in the stability analysis, the 
plateau instabilities vanish and only the ridge instabilities re- 
main. 

The 'ridge' instabilities with typical growth times of sev- 
eral thousand days owe their existence to a resonant coupling 
of wave modes with different wave numbers. As an example, 
the top panel of Fig. |S]shows the growth times of both Parker- 
type and background instabilities along a cut at i = 1 .64° in 




9.0 9.2 9.4 9.6 9.8 10.0 

B [1 0^ G] 



Fig. 6. Growth times r ( top) of Parker-type, ridge- and plateau- 
instabilities at A - 1.64° (see Fig. |5}- The ridge instabilities 
(grey shading) are due to resonant interactions of adjacent wave 
modes with m - 1 and 3 (middle) and m - and 2 (bottom), 
respectively, when the frequencies %((jj) are similar Dotted 
lines show the frequencies of the uncoupled modes in the corre- 
sponding single-star problem. Dashed and dotted-dashed lines 
indicate unstable modes with 3(to) + 0. 



Fig. |5] The panels below show the frequencies ''X(w) of the 
wave modes involved. In the single-star problem (dotted lines, 
dashed-dotted if unstable) the wave modes corresponding to 
different values of the azimuthal wave number, m, are indepen- 
dent. But the eigenmodes in the binary problem consist of a 
spectrum of coupled wave modes (see Fig.|5} with wave num- 
bers n - m,m + 2,m + 4, . . .: modes of wave number m with 
frequency %(a)) stimulate excitations of the wave modes m + 2, 
which, supposed %(a)) is close to the eigenfrequency of the ad- 
jacent wave mode, results in a resonant amplification and even- 
tually in the onset of an instability (solid lines in Fig.|6l dashed 
if unstable). 

In the following, we dismiss instabilities with very long 
growth times t > 10"^ d (^ 27 years), since they are not ex- 
pected to play any role for the eruption of magnetic flux tubes 
in active stars. 
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Fig. 7. Azimuth 0max of the radial envelope maxima. Due to 
the azimuthal 7r-periodicity only the maximum in the interval 
< 0inax < TT is shown, i.e., between the points Vi and on the 
line of centres implying the second maximum at 0max + The 
dashed framed region is the instability domain of eigenmodes 
with a dominating wave number m - 2. 




B [10*G] 

Fig. 8. Relative peak-to-peak variation, A||, |, of the radial en- 
velope. For ridge instabilities at low latitudes, A||,.| is close to 
unity, and even in the instability island at A ~ 15 . . .40° and 
low field strengths the resonant wave mode coupling leads to 
considerable values of A\^r\- 



4.3. Eigenfunctions 

The eigenfunctions, ^, inherit the property of the underlying 
problem and exhibit a ;7r-periodicity in azimuthal direction, im- 
plying a 0-dependent probability that a loop penetrates to the 
superadiabatic part of the convection zone and rises rapidly to- 
ward the surface. The preferred longitudes, 0n,ax (and 0max 
of loop penetration are determined by the maxima of the ra- 
dial envelope, l^rlmax = l|r(0max)l- They are shown in Fig. 
as a function of the equilibrium latitude Ao and field strength 
Bo- Background instabilities with growth times exceeding 10"^ d 
have been neglected, leaving behind only a fraction of a ridge 
at low latitudes and a plateau at A > 60° with t > 3 ■ 10"'d. 
This high-latitude feature is caused by the slow loss of equilib- 
rium due to the internal flow described above, which exhibits 
only small displacements in radial direction (here, llrl/ll^l ~ 
C?(10"^)). In contrast to the case of stable eigenfunctions, for 
which the envelope maxima are located exactly at the 'symme- 
try points', Qi 2 or Vi ^, the orientation of ||rlmax of unstable 
eigenfunctions show a phase shift with respect to these points. 
For instabilities at equatorial and low latitudes dominated by 
the wave number m - I, the maxima are located in the vicin- 
ity of Vi^2, i e., around the line connecting both stellar centres, 
whereas for modes with m - 2 they are distributed in a broad 
interval around the quadrature points 2i,2; preferred longitudes 
of loop penetration thus depend on the equilibrium parameters 
(Bq, Aq) as well as on the dominating wave mode. 

The radial envelope, is used as a measure for the prob- 
ability of loop penetration into the convection zone, subject to 
the normalisation ||r|max ■- 1- The larger the relative peak-to- 
peak variation, 

A||,| '^'•'■""^ ~ , (25) 

l^rlmax 

the larger is the possibility for the transition of a loop from the 
overshoot region to the convection zone around 0max- Figure 



|S] shows A||,.| as a function of Aq and Bo- For B > 10^ G the 
peak-to-peak variation is < 5% at intermediate latitudes and up 
to about 20% near the equator. Typically A||, | decreases with 
increasing field strengths since the stability properties strongly 
depend on /3S oc S/B^. High field strengths reduce the effect 
of the azimuthal (5-variation, so that magnetic flux tubes with 
smaller field strengths are, in general, more susceptible to tidal 
eff'ects. Particularly for ridge instabilities and instability islands 
at low field strengths, A\^r\ is significant and can even reach 
values close to unity, so that radial displacements are almost 
suppressed near the envelope minima. The approximation used 
in Eq. (I22> is not applicable here, since the dominance of the 
wave mode fr.o is not always given. Ridge instabilities are due 
to resonant interactions of adjacent wave modes (see Sec.|^} 
with comparable amplitudes while the amplitudes of all other 
modes are small. The case ||,.,o|, ||/-,2l s> VA: 0, 2, for ex- 
ample, yields the radial envelope 

ti" ex 1 + 2 ^^''f^f COS (20 + a2- ao) , (26) 

withao/a = ^"g|^r,o/2- For||r,ol ~ ||,-,2 1, the radial envelope drops 
to very small values at longitudes which are determined by the 
phase factors ao and a2- The low magnetic field strengths cor- 
responding to these particular instabilities imply, however, long 
growth times which are probably not relevant for the evolution 
of rising flux loops in stars. 

4.4. Parameter dependence 

We examine the dependence of the stability properties on tidal 
eff'ects in various binary systems which are characterised by the 
orbital period, T, the mass ratio, q, and the mass of the primary 
star. Mi, . Equations Q and (|5} show that the magnitude of the 
azimuthal variation caused by tidal effects is governed by the 
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Fig. 9. Comparison of growth times, t, for systems with var- 
ious orbital periods, T. The isoHnes mark the growth times 
T = 1,5, 10,50, 100 (thick line), 500, 1000 (dashed line), and 
5000 d. Instabihties with t > 10"* d are omitted. 



term 



47r2 



<1 



^ 10 



G 1+q M^T^ 

l+ql\RQl [MqI \d) 



-2 



(27) 



Here, the solar-type model is retained, which fixes and tq. 
The variation of Eq. J27I I due to a different value of ro is negli- 
gible since the equilibrium radius is restricted to the rather thin 
overshoot layer, which is prescribed by the stellar model. 

First we consider binary systems with orbital periods in 
the range T - 0.75 ... 10 d, mass ratio q - 1, and, follow- 
ing Eq. Q, separations a ^ 4.4 ... 25 Rq. Since the expansion 
parameter covers the range = 4.5 ■ lO""* . . .2.5 ■ 10"^, the 
analytical approximation for the stellar stratification, Eq. (|5j, is 
well applicable even for the shortest period. A comparison of 
growth times for different T is shown in Fig.|5] A smaller or- 
bital period has a stabilising effect, which is indicated by the 
shifting of corresponding isolines of growth times to higher 
field strengths. This effect is reinforced by the rotational flat- 
tening of the star since flux rings at a given, constant equilib- 
rium depth, ro, are embedded in slightly deeper and more sub- 
adiabatic layers of the overshoot region. The tidal effects on 
T are nevertheless small and primarily affect the background 
instabilities, leading to shifts of ridges and plateaus. With de- 
creasing system period, the high-latitude plateau grows toward 
smaller field strengths since the azimuthal variation of the den- 
sity and tidal force along the flux ring become larger; however, 
the growth times remain very large. The orientations of the ra- 
dial envelope maxima do not show a significant dependence on 
orbital period and are well represented by the reference case 
with r = 2 d (Fig. 0. The relative peak-to-peak variations, 
A\^r\, are shown in Fig. ^| A||r| increases significantly with 
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Fig. 10. Relative peak-to-peak variation, A\^r\, of the radial en- 
velope for different system periods, T. 

decreasing orbital period, which results in strongly preferred 
longitudes of loop penetrations to the convection zone in short- 
period systems. The tidal effects are very effective for flux rings 
with small field strengths, particularly for those which are li- 
able to ridge instabilities with A||, | close to unity. For Parker- 
type instabilities with short growth times, the relative peak-to- 
peak variation remains much smaller. Figure[TO]also shows that 
the peak-to-peak variation strongly decreases with increasing 
system period for T > 3 d and eventually becomes insignifi- 
cant. 

Considering a system with period T - 2d, we find that 
a variation of the mass ratio in the range q - 0.1 ... 10, has 
only marginal influence on the stability properties. In the do- 
main of Parker-type instabilities, the relative changes of the 
growth time and of the peak-to-peak variation with respect to 
the reference case are well below 0.1% and 10%, respectively. 
Furthermore, the orientations of the envelope maxima exhibit 
hardly any dependence on the mass ratio. 

In summary, the growth times show only a weak depen- 
dence on the binary parameters, except for the instability back- 
ground and the rotational stabilisation of flux tubes at short sys- 
tem periods. The latter effect is mainly due to the conservation 
of angular momentum and not innate to the tidal effects. The 
relative peak-to-peak variation of the radial envelope exhibits a 
considerable dependence on the orbital period T, while the de- 
pendence on the mass ratio q is almost negligible. Finally, the 
orientations of the envelope maxima are largely unaffected by 
variations of the binary parameters. 

5. Discussion 

Tidal effects and high rotation rates in close binaries cause 
characteristic changes of both the equilibrium and the stability 
properties of magnetic flux tubes stored in the convective over- 
shoot region. Owing to tidal effects, a flux ring in stationary 
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equilibrium experiences periodic variations of its environment 
in azimuthal direction, which modify its characteristic eigen- 
modes. To first order, the azimuthal variation of a given quan- 
tity depends on its local scale height. Since the stability proper- 
ties are dominated by the superadiabaticity, which is the most 
strongly depth-dependent quantity in the overshoot region, tidal 
effects are mediated mainly by the azimuthal variation of 6. 

The growth times of Parker-type instabilities with t < 
1000 d are only marginally affected by the presence of the com- 
panion star, whereas for instabilities with t > lOOOd the influ- 
ence of tidal effects is considerably stronger. This applies par- 
ticularly to the instability background, since ridge and plateau 
instabilities owe their existence to the binary character and are 
most susceptible to changes of the tidal effects. However, since 
the eigenfrequency is a global parameter describing the whole 
flux tube and thus representing an average over the properties 
of all tube segments, the dependence of co on the binary pa- 
rameters is dominated by the axially symmetric influence of 
stellar rotation and not so much by the non-axially symmetric 
tidal interactions: for axially symmetric changes of the equi- 
librium configuration, e.g., those due to the rotational flatten- 
ing, the changes are identical for each tube segment and thus 
may add up to a considerable deviation of oj, whereas for non- 
axial, periodic changes along the flux ring, like those due to 
tidal effects, the changes of individual tube segments mutually 
cancel and eventually average to only minor deviations of w. 
The growth times of background instabilities are typically sev- 
eral thousand days and more and therefore very long compared 
to the system period as well as to the growth times which are 
relevant for rising flux tubes in the Sun. Whether these insta- 
bilities are important for loop formation at preferred longitudes 
seems questionable, although they possess the most conspicu- 
ously asymmetric eigenfunctions. 

In contrast to the modest effect on the eigenfrequencies and 
growth times, tidal effects can significantly alter the azimuthal 
structure of the eigenfunctions since the periodic variations 
along the tube affect the displacement of each segment indi- 
vidually. The resulting eigenfunction, ^(<p), characterised by its 
radial envelope ||, |, consists of a superposition of coupled wave 
modes and inherits the 7r-periodicity of the azimuthal variation 
of the equilibrium properties. In contrast to stable eigenmodes, 
whose envelope maxima are located either at the line of centres 
or at the quadrature points, the maxima of unstable eigenmodes 
typically exhibit a phase offset with respect to these points; this 
result is schematically illustrated in Fig. ^2 The orientations 
of the maxima are not unique but depend on the equilibrium 
configuration, predominantly on its latitude and on the wave 
number of the dominating wave mode. Considering Fig.|2l one 
might expect that the quadrature points would be favoured for 
rapid loop growth, since they are located in less stable layers 
of the overshoot region and nearest to the superadiabatic con- 
vection zone. However, the Parker-type instability is driven by 
a downflow from the loop summit, which increases the density 
contrast with respect to the environment. Following the equa- 
tion of continuity, a net downflow implies an asymmetry of 
mass flux at the loop summit, i.e., the amount of plasma flowing 
down in the preceding slope of the loop outbalances the mass 
flux coming up in the following loop slope. The azimuthal vari- 




Fig. 11. Schematic illustration of flux tubes in a single star 
{left) and in a binary star (middle and right). For the latter, the 
direction to the companion star is indicated by the central ar- 
row. The dashed line indicates the equilibrium flux ring, the 
grey line an actual displacement of the perturbed flux tube (here 
for a wave mode with m =1), and the solid black lines the enve- 
lope, llrl, of the eigenfunction. In the middle panels, two stable 
modes are shown with the maxima of their radial envelope lo- 
cated either at the quadrature points [upper panel) or at the line 
of centres [lower panel); unstable modes in binaries typically 
exhibit phase offsets with respect to these points [right). 



ation of the internal equilibrium density and velocity, Eqs. 
and (II 1> . enhance the asymmetry in both slopes and thereby 
favour maximal displacements at intermediate longitudes. In 
contrast, the extrema of stable eigenmodes prefer either the 
quadrature points or the directions along the line connecting 
both stellar centres since at these 'symmetry points' (concern- 
ing the internal equilibrium density and velocity structure) the 
downflow and thus buoyancy effects can be minimised. 

The peak-to-peak variation, A\^r\, serves as a measure for 
the probability of loop penetration into the convection zone, 
based on the assumption that the loop segments lifted furthest 
toward the superadiabatic region are the most likely to leave 
the overshoot region. The radial envelope determines the dis- 
placement of the perturbation wave and thus the formation of 
rapidly growing loops inside the convection zone is favoured 
around its maxima at <p,^.^^ and 0n,ax + ^- In the case of close 
binaries with orbital periods of a few days, the peak-to-peak 
variation of fast growing Parker-type instabilities is of the or- 
der of several percent, up to about 20%. In the case of ridge 
instabilities, radial displacements are strongly suppressed at 
some longitudes, giving rise to distinct preferred longitudes 
in the vicinity of ^^ax and 0n,ax + ^- Although it can not be 
excluded that ridge instabilities contribute to the non-uniform 
spot distribution, their importance for the formation of pre- 
ferred longitudes at the stellar surface is questionable since 
the relevant equilibrium configurations are restricted to rather 
small parameter domains (Z?o, Aq), where instabilities have 
large growth times. Based on the parameter study in Sect. 14.41 
we would predict that preferred longitudes appear (in bina- 
ries with solar-type components) for orbital periods J < 5 d. 
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which is in agreement with the results of iHeckert & Ordwavl 
( I1995I) . However, in some binaries with longer periods (consist- 
ing of gia nt components) preferr ed longitudes have also been 
observed dBerdvugina & TuominerLl 998). 

The linear stability analysis establishes the existence of 
preferred longitudes at the bottom of the convection zone. 
However, the determination of the actual longitude where an 
individual loop enters the convection zone is not possible since 
our procedure only allows statistical predictions. The tendency 
to constitute preferred longitudes is instead discernible for a 
large ensemble of similar flux tubes, supposed that the environ- 
ment conditions in the overshoot region, e.g., the amount and 
distribution of overshooting convective motions, are uniform 
and unaltered for a sufficiently long time. Preferred longitudes 
at the bottom of the convection zone are not necessarily sim- 
ply related to the observed preferred longitudes of starspots at 
stellar surfaces, because the evolution of rising loops through 
the convection zone is still subject to tidal effects^. The pre- 
ferred longitudes of flux eruption at the stellar surface based on 
the flux tube model described above have to be determined by 
non-linear numerical calculations. Results of such calculations 
are presented in the subsequent paper II. 

6. Conclusion 

The equilibrium and stability properties of toroidal magnetic 
flux tubes in the overshoot region of an active component of a 
close binary system are considerably affected by tidal effects. 
The breaking of the axial symmetry due to the presence of the 
companion star leads to azimuthal variations along the tube of 
the equilibrium configuration and of the probability of the pen- 
etration of unstable loops to the convection zone. This demon- 
strates the existence of preferred longitudes at the bottom of the 
convection zone, which are well pronounced for orbital periods 
of a few days but strongly fade with increasing period. The ori- 
entations of the preferred longitudes depend on the equilibrium 
latitude, the field strength, and the dominating wave mode of 
the unstable eigenmode. Despite their rather small magnitude, 
tidal effects are nevertheless capable of influencing the dynam- 
ics of magnetic flux tubes considerably because they are en- 
hanced by the sensitive dependence of the stability properties 
on the superadiabaticity of the local environment and by the 
resonant coupling between interacting wave modes. 
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Appendix A: Numerical treatment 

The perturbations of all physical and geometrical quantities 
(magnetic field, flow velocity, density, position, curvature,. . . ) 
arising from the small perturbations of an equilibrium flux ring 



^ If the tube evolution would be linear throughout the convection 
zone, the maxima of the radial envelope would, in fact, also determine 
preferred longitudes of flux eruption at the surface. 



are expressed in terms of the displacement vector ^{(p) (see, 
e.g.,|Ferriz-Mas & Schiissler 1993). The substitution of the cor- 
responding expressions in the equation of motion followed by a 
linearision with respect to ^ yields Eq. J14l i. with the coefficient 
matrices^ 

M„ = 



1 1 1 



Ml 



4Q2 Ml 

ri^-iW + ri^z 



1 1 1 -H 
J. i J. -r ^^^a_^ 



2 n Mi 
1 



(A.l) 
(A.2) 

(A.3) 

(A.4) 



8t 2gn gb 

gb } 



8* 



gtgt glgn gtgb 
gngt gngn gngb 
gbgt gbgn gbgb 



gn gb 


j 



where the constant matrices 
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0] 




' 





1 ^ 







1 







-1 





















lo 





ij 




[ 





oj 




l-l 





oj 



and the abbreviations 
1 Pe 1 



A 



(1 - Ml) p^y K 

-r^[{\-V)Ml + [\-Ml)v\y^ 

s* 

2 

P6- 



rf^A(l-M^) 



(A.5) 

,(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A. 10) 



have been used. The components - {guS ■ x) ,x e{t, n, b] of 
the effective gravitational acceleration geff are expressed in the 
co-moving trihedron. 

Owing to the azimuthal variation of the equilibrium quan- 
tities induced by the presence of the companion star, the coef- 
ficient matrices iA.H - iA.5\ can be expressed in the form of 
Eq. il5i and comprise small contributions which are, in low- 
est order, 7r-periodic in the longitude <p. Following the ansatz in 
Eq. ilH . these contributions cause the coupling of wave modes 
In with wave numbers n = m, m ± 2, m + 4, . . ., which is ex- 
pressed either by the 3-term recursion formula in Eq. ( I18> or 
by the Hill-type algebraic system 



■Cm-IK C 



-IK 



o 



-Cm 



Cm 



lo 



o 



m+2K C„,+2K 



•+2K ) 



= 0, (A.ll) 



In the case /3 » 1 . 
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with K oo and the block matrices 

-cJ{M„,c-iM„^s)\, (A. 12) 

K = ^ [(« + 2) (/Al^.c - Al^,.5 - w (A1^,,c + 

-a? {M„,c + iMn^s)] , (A. 13) 

and 

C„ - -n^& - nojM^ifl - ci/M„,q + M^,o 

(nM^.o + wM,o) ■ (A. 14) 

The eigenfrequencies oj are determined by the roots of the 
dispersion relation, det ['7Y(w)] = 0. A treatment of the case 
K oo is not practical, but since the coupling between ad- 
jacent wave modes is small [X.,'R ~ (9(e^)], the wave mode 
spectrum can be truncated at a finite number of constituents, 
here A' = 3. A verification of the eigenfrequency and the deter- 
mination of the eigenfunctions ^{(p) is based on the following 
approach. Let the independent solutions x^i^{<p), i - 1, . . . , « of 
a system 

x'^^(^)x with ^(^ + 7t) ^ Ji((p) (A. 15) 

represent the fundamental matrix X{<p) - . . . with 
det[<Y(0)] 7t 0. Furthermore, we have X((f) + In) = -Y((^)A1, 
where M. is the constant transition matrix with det(AI) + 0. 
Since an arbitrary solution of Eq. jA.15> can be represented by 
a constant vector c in the form x(<p) = P((<p)c, it follows that 
x(<p + 2n) - X{(p + 2n)c - X{(p)Mc. Periodic solutions with 
x{4> + 27r) = x{^) thus require Mc - c and an eigenvalue of the 
transition matrix M equal to unity. If the fundamental matrix 
fulfils the initial condition X{Q) - £, the transition matrix is 
calculated by integrating over a period of 2jt, i.e., X(2jt) - M. 
In our case, the transition matrix M and its eigenvalues depend 
on the eigenfrequency, cj. If an eigenfrequency determined by 
the dispersion relation above does not entail Mc - c, w is 
iteratively improved and the eigenfunction ^{^) determined. 
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